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Using our new photometric and spectroscopic observations as well as all available published data, we present a new interpretation 
of the properties of the peculiar emission-line binary UX Mon. We conclude that this binary is in a rare phase of fast mass transfer 
between the binaiy components prior to the mass ratio reversal. We firmly establish that the orbital period is secularly decreasing at a 
rate of P = -0.260 seconds per year. From several lines of reasoning, we show that the mass ratio of the component losing mass to 
the mass-gaining component q must be larger than 1 and find our most probable value to be g = 1.15 ± 0.1. The BINSYN suite of 
programs and the steepest descent method were used to perform the final modeling. We modeled the star as a W Ser star with a thick 
disk around its primary. Although the remaining uncertainties in some of the basic physical elements describing the system in our 
model are not negligible, the model is in fair agreement with available observations. Only the nature of the light variations outside the 
primary eclipse remains unexplained. 
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1. Introduction 

UX Monocerotis (HD 65607, SAO 135333) is an eclipsing bi- 
nary with an orbital period of 5'!90. Its binary nature was first 
reported by Woods & Shapley (1928), who also noted that its 
brightness at maximum light is not constant. Gaposchkin (1947) 
studied 32 spectrograms from the McDonald Observatory and 
1389 photographic plates from the Harvard Observatory. He 
identified the star eclipsed at the deeper minimum (hereafter, 
the primary) as an A7p object and the other component (here- 
after, the secondary) as a G2p star He arrived at surprisingly 
low masses of Mi = 0.73 Mq and M2 - 0.74 Mq for the pri- 
mary and the secondary, respectively. Struve (1947) analyzed 
152 photographic spectra secured in 1947. He obtained radial 
velocity (RV hereafter) semi-amplitudes of Ki - 140km s"' and 
K2 = 60 km s ' for the primary and secondary, respectively. This 
implies a more massive secondary, M2 sin"' / = 3.4 Mq, than the 
primary. Mi sin^* / = 1 .5 Mq, which would mean that UX Mon is 
in the rare evolutionary stage prior to the mass reversal (cf . e.g. 
Crawford 1955; Morton 1960). Struve also discussed the ori- 
gin of the Balmer emission, which he identified as circumstellar 
matter located between the two stars. 

Hiltner et al. (1950) published 486 yellow photoelectric ob- 
servations with an effective wavelength of 5300 A. They con- 
firmed the light variations outside the primary eclipse. They also 
noted that the durations of both eclipses are the same. Wood 
(1957) obtained light curves of the system in the yellow and blue 
filters. He observed short (hours) and long (six or more days) 
non-orbital light changes. Lynds (1956, 1957a,b) secured and 
analyzed 333 UBV observations. He derived the ratio of radii, 
k = 0.522, the radius of the secondary relative to the orbital sep- 
aration, r2 - 0.358, and the inclination, / = 83°6. 
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Scaltriti (1973, 1976) concluded that the primary is a 5 Scuti 
variable and that there had been an abrupt change of the orbital 
period in the past. 

Umana et al. (1991) found that the UX Mon system is a 
microwave source and concluded that the microwave radiation 
originates from the gyrosynchrotron radiation of electrons in a 
magnetic field. 

Olson & Etzel (1995) carried out a detailed photometric and 
spectroscopic study of UX Mon based on their new ubvyl ob- 
servations and a series of electronic spectra. They derived the 
most accurate RV semi-amplitude of the secondary yet obtained, 
K2 =108.3+1.9 km s"\ but were unable to detect any measur- 
able spectral lines of the primary. In contrast to Struve (1947), 
they inferred a mass ratio secondary /primary of <7 = 0.8, based 
on light curves analysis. Their model also led to the surface tem- 
peratures Ti a: 8000 K and T2 ~ 5500 K for the primary and 
secondary, respectively. They concluded that the gas stream from 
the secondary hits the primary directly, without forming an ac- 
cretion disk around it. 

Analyzing published photometry and spectroscopy, Ondfich 
& Harmanec (2003) concluded that the orbital period of 
UX Mon has been secularly decreasing and provided new 
quadratic ephemeiis. They also discussed the emission lines of 
highly ionized atoms in the HST/GHRS spectra indicating that 
UX Mon is a W Ser star W Ser stars are strongly interacting bi- 
naries characterized by the presence of numerous emission lines 
in their UV spectra, which are believed to originate in an ex- 
tended gaseous envelope around their mass-gaining components 
(cf. e.g. Plavec & Koch 1978; Plavec 1980, 1992). However, 
the mass transfer rate (10^** - 10"^ M© y"'), which Ondfich 
& Harmanec (2003) derived for UX Mon from the quadratic 
ephemeris, appears to be too low compared to the typical val- 
ues of 10"^ M0 y"' measured for W Ser binaries (Plavec 1980, 
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Table 1. Journal of photometric observations of UX Mon 



Sources 


Time interval 


No. of. 


Passbands 




(HJD-2400000) 


obs 




(1) 


18600.0-18605.7 


38 


VIS 


(2) 


33246.9-33381.6 


486 


^5300 


(3) 


33300.8-33399.7 


401 


i5200, ^4250 


(4) 


35084.9-35153.9 


999 


UBV 


(5) 


41626.7-41764.3 


875 


i5150, V 


(6) 


47181.7-48656.9 


5056 


uvbyl 


(7) 


47964.6-48938.2 


76 


V 


(8) 


51868.8-53896.5 


268 


V 


(9) 


52655.4-54219.3 


402 


UBV 


(10) 


52655.4-54219.3 


137 


vis 


(11) 


51448.0-54228.6 


1453 


BV 



References. (1) Woods & Shapley (1928); (2) Hiltner et al. (1950); (3) 
Wood (1957); (4) Lynds (1957a); (5) Scaltriti (1973); (6) Olson & Etzel 
(1995); (7) Ferryman & ESA (1997); (8) Pojmadski 1997; (9) Hvar 
observations (this paper); (10) Meyer (2007); (11) Olson et al. (2009). 



1992). Independently, the secular period decrease was also noted 
by Meyer (2006). 

Olson et al. (2009) published a new study of UX Mon based 
on systematic photometry secured over a period of eight years. 
They found almost the same system parameters as in their previ- 
ous paper (Olson & Etzel 1995). Focusing mainly on light vari- 
ations outside the primary ecUpse, they concluded that these are 
probably caused by the variable mass transfer. 

2. Observations 

2.1. Photometry 

Altogether, 22 dilFerent datasets of photometric observations of 
UX Mon from 1 1 different sources were used. Table 1 lists basic 
observational data. More information about the datasets used in 
this paper is given below. 

(1) Photographic photometry, consisting of 38 phase-averaged 
normal points published by Woods & Shapley (1928) only 
in orbital phases. The accurate dates of observations are un- 
known limiting their usefulness in our study. 

(2) These observations were obtained at McDonald Observatory 
from November 26, 1949 to April 10, 1950. 

(3) Wood (1957) obtained 401 photoelectric observations, di- 
vided into two {green and blue) datasets. 

(4) There are 333 individual observations in each of the UBV 
passbands in this dataset, obtained during twenty nights be- 
tween December 7, 1954, and February 15, 1955. The instru- 
mental color system is essentially the same as the standard 
UBV system, except for the zero point. 

(5) Two datasets of green and yellow photoelectric observations 
contain a total of 875 individual observations. 

(6) These five sets of uvbyI{Krori) observations consist of 5056 
observations in total. They were secured between January 
21, 1988, and February 4, 1992 and analyzed by Olson & 
Etzel (1995). Upon our request. Dr. Olson kindly allowed us 
to have copies of these observations at our disposal. These 
data were subsequently published by Olson et al. (2009). 

(7) The Hipparcos Hp photometry of UX Monocerotis was ex- 
tracted from the data archive published by Ferryman & ESA 

(1997) and transformed to the standard Johnson V magni- 
tude using the transformation formulae derived by Harmanec 

(1998) and Harmanec & Bozic (2001), which are based on 



numerous all-sky UBV observations of many stars observed 
at Hvar. 

(8) All-Sky Automated Survey (ASAS hereafter) is an observa- 
tional program for monitoring stars of the southern hemi- 
sphere brighter than x!l4th magnitude (Pojmaiiski 1997, 
2001). Altogether, 268 Johnson V observations of UX Mon 
were secured between November 20, 2000 and June 10, 
2006. 

(9) Altogether, 402 new photoelectric UBV observations were 
secured during two seasons with the 0.65-m Cassegrain re- 
flector of the Hvar Observatory. The transformation from 
the instrumental to the standard Johnson system was carried 
out by the program HEC22 (Harmanec et al. 1994), which 
uses non-linear transformation formulae. The most devel- 
oped version rel.16.2 was used, which allows modeling of 
variable extinction during the observing nights. 

(10) This dataset consists of visual estimates by Meyer (2006, 
2007). 

(11) Olson et al. (2009) obtained 1453 individual measurements 
in Johnson's B and V filters with the 0.4m APT at Fairbom 
Observatory in southern Arizona, USA. 

2.2. Spectroscopy 

In total, 39 new electronic spectra obtained at three observatories 
are included in this paper. 

(1) Six spectra were obtained at the coude focus of the Ondfejov 
Observatory 2.0-m reflector with a linear dispersion of 17 A 
mm"^ (two-pixel resolution 12600), which cover a spec- 
tral range 6300-6700 A. Subtraction of bias, flatfielding, 
creation of 1-D spectra, and wavelength calibration has 
been routinely carried out shortly after the observations by 
Dr M. Slechta. Rectification, heliocentric correction of time 
and radial velocity, and the removal of cosmic spikes was 
carried out by DS with the help of the program SPEFO (Hom 
et al. 1996; Skoda 1996) in its form developed by the late 
Mr. J. Krpata. 

(2) Thirteen spectra were obtained in the coude focus of the 
1 .22-m reflector of the Dominion Astrophysical Observatory 
(DAO hereafter). Their linear dispersion is 10 A mm ' (two- 
pixel resolution 21700), while the range of wavelengths is 
6100-6700 A. The initial reduction of spectra (bias removal, 
flatfielding, and creation of 1-D spectra) was carried out by 
SY using IRAF. The wavelength calibration and all subse- 
quent reductions were carried out by DS in SPEFO. 

(3) Twenty echelle spectra were obtained at the coude focus of 
the 2.0-m reflector of TLS Tautenburg. These spectra have 
the best hnear dispersion of all three sets: 3.2 A mm ' (two- 
pixel resolution 63000). They cover the spectral range 4700- 
7085 A. All reductions, including the rectification and he- 
liocentric corrections were carried out by HL using stan- 
dard MIDAS packages and a special routine to determine the 
nightly instrumental RV zero points from a large number of 
telluric O2 lines in the spectra. 

3. Analysis 

3.1. Secular changes of the orbital period 

Ondnch & Harmanec (2003) proposed that there had been a sec- 
ular decrease of the orbital period. We decided to re-investigate 
their finding, using more photometric observations than they had 
at their disposal. 
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Table 2. Times of tlie primary minima and values of the instantaneous 
orbital period as derived by FOTEL assuming locally linear ephemeris. 
The column O-C2 gives the differences with respect to quadratic 
ephemeris (Eq. 4). 



Source 


Epoch 


Period 


O-C2 




(HJD-2400000) 


(days) 


(days) 


(2) 


33246.1720(092) 


5.905882(563) 


-0.0053 


(3) 


33299.3269(063) 


5.904718(520) 


+0.0082 


(4) 


35088.3882(047) 


5.904210(689) 


-0.0208 


(5) 


41630.6478(157) 


5.904829(920) 


-0.0074 


(6) 


47180.8866(007) 


5.904494(004) 


-0.0028 


(8) 


51869.0342(078) 


5.904404(038) 


+0.0016 


(9) 


53746.6472(038) 


5.904432(067) 


+0.0042 


(11) 


52772.4039(019) 


5.904378(013) 


-0.0008 



Table 3. The quadratic ephemeris, based on all photometric datasets of 
Table 2, derived with FOTEL. 



Quantity 


Value 


7^prim.mm. [HJD] 


2452666.1319 ±0.0012 


Po [days] 


5.9044365 ± 0.0000020 


f'T=To [days per day] 


(-8.23 ± 0.18) xlO-» 


a [days] 


(-2.43 ± 0.05) X 10-8 



We first calculated the local linear ephemeris of the sources 
listed in Table 1 with a sufficient number of measurements by 
using the latest publicly available version of the program FOTEL 
(Hadrava 2004a). We allowed for the convergence of the epoch 
of the primary minimum and the local period. The results are 
summarized in Table 2. We also added a third column, O-C2, 
listing differences between the local epochs and those calculated 
from the quadratic ephemeris (see below). 

In FOTEL, there is also an option that allows us to de- 
rive a period change, P, that exactly corresponds to a quadratic 
ephemeris, as one of the unknowns of the solution. As discussed 
by Harmanec & Scholz (1993), a quadratic ephemeris given by 



T = To + Po-E + a-E^ 



(1) 



holds for the times of minima, where Po denotes the value of the 
period at the time of reference minimum Tq. Treating the epoch 
£ as a real number characterizing the cycle and phase, and de- 
noting T as the time corresponding to E (i.e. to the given orbital 
phase at that cycle), one can write the instantaneous orbital pe- 
riod, P, as 



dT 

P= — ^Po + 2aE. 
dE 

The instantaneous change of the period, P, is then 

. dE 2a 
P = 2a— = — . 
dT P 



(2) 



(3) 



This shows that the quadratic ephemeris leads to a Unear change 
of the period with epoch but not with time. For instance, for 
negative values of a, the rate of the period change increases with 
time. 

We therefore used the corresponding option in FOTEL to 
derive the quadratic ephemeris given in Table 3 by combining 
all datasets from Table 2. 

The explicit expression for the quadratic ephemeris is 




rmin.i = HJD2452666.1319 + 5.9044365E 



■ 2.43 X IQ-^E^, 



(4) 



-3000 -2000 
epoch 

Fig. 1. 0-C in days with respect to the linear ephemeris. Crosses with 
error bars denote the times of minima from Table 2. The empty squares 
denote epochs of minima of sources (1) and (10). The times of individ- 
ual observations for source (1) are unknown, and measurements from 
source (10) are of insufficient acciu-acy. Quadratic ephemeris is shown 
with a solid curve. 



where the coefficient of the quadratic term is calculated as 
PqPt=ToI'2- The change of the orbital period corresponds to 
Pt=To - -0.260 seconds per year. 

The local epochs of minima from Table 2 were compared to 
the global linear ephemeris and the corresponding 0—C devia- 
tions are shown in Fig. 1 . They clearly show the parabolic change 
indicating a secular decrease of the orbital period, in accordance 
with the findings of Ondfich & Harmanec (2003). 

Note that the epoch of the early photographic observations 
by Woods & Shapley (1928) does not fit well with the extrap- 
olated prediction of the quadratic ephemeris (Eq. 4), but it also 
confirms that the period has been decreasing. Since the dates of 
these early observations are not known, we could not include 
this dataset into our FOTEL solution to obtain yet more accu- 
rate values for the quadratic ephemeris and the rate of the period 
change. 

In Fig. 2 we fold all photometry near the primary minima 
from the datasets of Table 2 with the quadratic ephemeris (Eq. 4), 
obtaining a very coherent fit. Although the depths of the minima 
differ, of course, for each passband, the phase of the primary 
mid-eclipse is the same for all of them. 

Note that Ondnch & Harmanec (2003) derived a smaller 
value of the period decrease, P = (-2.49 + 0.18) x 10 days 
per day. For a fully conservative mass transfer, they estimated 
M = 1.4 X IQ-^ Mq per year. We note, however, that their value 
refers to a mass transfer rate per one day, and that the correct 
mass transfer rate per one year should be 0.52 /zMoy ' , i.e. typi- 
cal of strongly interacting binaries (Plavec 1980). Our new result 
imphes a mass transfer rate that is even higher than the value by 
Ondnch & Harmanec (2003) (see below). 



3.2. Spectroscopic mass ratio 

The mass ratio of the binary systems with a Roche-lobe fill- 
ing component can be estimated from the (plausible) assump- 
tion that the rotation of this star is synchronized with the orbital 
revolution (Andersen et al. 1989; Harmanec 1990). The essence 
of the method lies in the relative dimensions of the Roche lobe 
depending solely on the binary mass ratio q, while the absolute 
radius of the spin-orbit synchronized Roche-lobe filling star is 
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Fig. 2. A phase plot of all photometric observations around the phase of 
primary minimum using quadratic ephemeris (Eq. 4). 



Fig. 3. )^ value as a function of V2 sin i for various values of f,. A 95% 
confidence level is shown with a dashed line. 



given by the equatorial rotational velocity, the incUnation of the 
rotational axis (assumed to be the same as the orbital inclination) 
and by the rotational period, which is identical to the binary or- 
bital period. The third Kepler law is used to derive the binary 
separation. 

To apply a similar procedure to UX Mon, we measured the 
rotational velocity of the secondary, V2 sin using a mean spec- 
trum constructed from three spectra closest to the phase of the 
primary mid-echpse (phase ^ = 0) and compared it with syn- 
thetic spectra. For that, a grid of model atmospheres was con- 
structed using the LLmodels2 program (Shulyak et al. 2004), 
which allows us to calculate 1-D plane-parallel, LTE model at- 
mospheres for stars with effective temperatures of 5000 K and 
higher. Line absorption is taken into account directly for each 
line taken from a line list selected by the user. For the spectrum 
synthesis, we used the program SyiJthV (Tsymbal 1996). 

To check whether the assumption of spin-orbit synchronism 
is in agreement with the value ofq = 0.8 obtained from photom- 
etry alone (Olson & Etzel 1995; Olson et al. 2009), we adopted 
the values derived by the same authors: the temperature of the 
secondary T2 - 5500 K and its surface gravity acceleration log g 
= 2.9. Free parameters of the fit were microturbulent velocity, 

and V2 sin i. We also allowed for the adjustment of the lo- 
cal continuum of the observed spectrum to that of the synthetic 
one, compensating for possible large-scale misalignments due 
to an imperfect continuum placement during the spectrum re- 
duction. To avoid the broad Balmer and stronger telluric lines, 
the analysis was restricted to the 4897 to 5777 A wavelength 
range. The best fit was obtained for - (0.68 + 0.23) km s"' 
and V2 sin i = 90.7 ±3.0 km s The^^ values as a function of 
V2 sin i for various values of ^1 are shown in Fig. 3. 

We used the value of V2 sin/ - 90.7 + 3.0 km s"^ and the 
orbital period P = 5.9044365 days for the reference epoch de- 
rived in the previous section (first term in Eq. 4). In addition, 
we used the inclination i = 89.5^ and the RV semi-amphtude of 
the orbital velocity of the secondary K2 = 108.3 km s"' derived 
by Olson & Etzel (1995). Following the numerical procedure 
presented in the appendix of Harmanec (1990) (see also the be- 
ginning of this section), we arrived at ^ = 1.3. This estimate of q 
difl'ers substantially from the value of ^ = 0.8 derived by Olson 
& Etzel (1995). Note that we estimate the mass ratio again with 
the same method but using new values of all relevant parameters 
derived in this paper at the beginning of Sect.4. 





































i\ 













0.4 0.6 
orbital phase 



Fig. 4. The RV curve of the secondary given by Eq. (5) is shown by 
a dashed sinusoidal curve. The individual mean RVs of the secondary, 
based on SPEFO measurements, with their error bars are also shown. 



To find a more accurate value of the mass ratio, q, we used 
all of the spectra at our disposal (cf. Sect. 2.2). The RVs of in- 
dividual unblended spectral lines were measured with SPEFO, 
which displays direct and reverse traces of the Une profiles su- 
perimposed on the computer screen and the user can slide them 
to achieve their exact overlapping for the studied detail of the 
profile. The orbital phases were calculated with the quadratic 
ephemeris (Eq. 4). Many suitable strong lines of the secondary 
could be identified. The RV curve of the secondary was calcu- 
lated by performing a non-linear least squares fit assuming a cir- 
cular orbit 



RV2 = (104.3 ± 1.1) sin(2;r0) km s'^ - 0.2 ± 0.8 km s" 



(5) 



The resulting theoretical RV curve is plotted with the individ- 
ual RV measurements in Fig. 4. A systematic deviation of the 
measured RVs from the calculated ones is visible before the sec- 
ondary eclipse (<p - 0.5). We attribute this difference to the 
Rossiter-McLaughlin effect. We return to this question in later 
sections. 

To find spectral Unes originating in the atmosphere of the 
primary proved to be a much more difficult task. A comparison 
of the spectra obtained near phases (pi - Q and 02 = 0.5 re- 
vealed that the most promising candidates are the following iron 
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600000 
595000 
590000 
585000 
580000 
575000 



Kg [km s" ] 



Fig. 5. a) Sum of squares of residuals, ssr, calculated with KOREL as 
a function of converged values of radial velocity semi-amplitude of the 
secondary, K2, and inverse mass ratio 1 /q shown as a map plot, b) The 
same as in a) but with l/q axis scaled for better clarity. 



lines: Fe II 4923.921 A, Fe II 5018.434 A, and Fe II 5169.030 A. 
However, the direct RV measurements with SPEFO, even for 
spectra from which the secondary lines with proper velocity shift 
were subtracted, resulted in a very uncertain RV curve display- 
ing a large scatter. We, therefore, attempted the spectral disen- 
tangling using the latest publicly available version of the pro- 
gram KOREL (Hadrava 2004b). We concentrated on finding the 
best values of two parameters: the RV semi-amplitude of the 
secondary, K2, and the mass ratio, q. We derived many solu- 
tions mapping the plausible range of both parameters. Since the 
value of K2 is fairly well constrained by the solution given in 
Eq. (5) and by an independent result by Olson & Etzel (1995) 
(K2 = 108.3 km s '), we only inspected the K2 in the range 
from 100 to 1 12 km s~' . On the other hand, we treated the mass 
ratio as unconstrained over a wide range of values between 0. 1 
and 10. We allowed convergence only for the time of the primary 
epoch and kept all other parameters fixed. 

In Fig. 5, we plotted the sum of squares of the residuals, 
ssr, as a function of the two parameters K2 and <7 as a map plot. 
Although the Fig. 5 represents only a crude mapping of the pa- 
rameter space, it is quite clear (see Fig. 5a)) that the best-fit so- 
lution lies in the region 0.5 < l/q < 1.0. 

To obtain more accurate values of the parameters K2 and q, 
we selected 20 initial parameter sets with the lowest value of 
ssr and allowed for the convergence of K2 and q while keeping 



Fig. 6. The disentangled spectra of the two components. 



everything else fixed. We obtained K2 = 104.6 km s^'and l/q = 
0.87 as the values of the lowest ssr. By estimating the error in 
1/g to be 0.05 and taking the error in K2 from Eq. 5, we end up 
with 

K2 = 104.6+1.1 kms 
q = 1.15 + 0.05, 
Ki = 120.6 ±5.4 km s"^ 



Ml sin^ i = 3.24 ± 0.38 Mq, 
M2sin^' = 3.74 + 0.28 Mq, 
flsin; = 26.29 + 0.64/^0. 



(6) 



Fig. 6 shows the disentangled spectra of both components. 
The spectrum of the primary is shifted upwards by 0.2 for clarity. 
The three strong Fe II lines are clearly visible in the spectrum of 
the primary. 

In the approximation of fully conservative mass transfer, a 
period decrease in time leads to the mass ratio q > I. Ondrich 
& Harmanec (2003) carried a similar analysis and came to the 
same conclusion, challenging the value of q - 0.8 published by 
Olson & Etzel (1995). Our new value of ^ = 1.15 agrees with 
their analysis. 

Using the well-known formula for fully conservative mass 
transfer (cf e.g. Harmanec & Scholz (1993)) and our estimates of 
the masses from the RV measurements, we obtain a mass transfer 
rate of M = (-3.6± 1.3) fiMgy'^. This value for the mass transfer 
rate is typical of W Ser stars (Plavec 1980, 1992). 

We also note that the region around Ha exhibits a strong 
emission at all phases (Fig. 7). All 39 available spectra are 
shown. One can see that changes in the Ha profile are linked 
to orbital motion of the primary and secondary, shown in Fig. 7 
by solid and dotted sinusoids, respectively. 

3.3. Light curve modeling 

Considering that mass transfer is present in the system and that 
Ondrich & Harmanec (2003) found the emission lines of Si IV 
around 1400 A and of C IV around 1550 A in the UX Mon spec- 
tra, we felt encouraged to model the system as a W Ser star. This 
required us to model an optically and geometrically thick disk 
surrounding the primary star. 

We used the BINSYN suite of programs (Linnell 1984; 
Linnell & Hubeny 1996; Linnell 2000) with a parameter opti- 
mization procedure based on the steepest descent method devel- 
oped by DS. The optimization procedure controlled the flow of 
BINSYN software, automating the process. 

Using the results from the previous section, we kept q - 1.15 
and Ml sin^ ; = 3.24 M© constant. As can be seen from the light 
curve, the eclipse is total or close to totaUty. Therefore, we can 
estimate the mass of the primary to be in the range of 3.24 Mq 
to 3.4 Mq. Assuming that the primary is a main-sequence star, 
we estimate its radius Ri ~3 Rq and temperature Ti ^ 13000 K 
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6550 6555 6560 6565 6570 6575 
wavelength 

Fig. 7. A plot of all 39 Ha profiles available to us, arranged according 
to their orbital phases (via appropriate shift of their continua along the 
ordinate labeled with the orbital phases). The solid and dotted sinusoids 
represent the orbital motion of the primary and secondary, respectively. 



(Harmanec 1988). Note that this assumption probably underesti- 
mates the actual radius because such large stars in close binaries 
are probably inflated with respect to their main-sequence coun- 
terparts. 

Although the effect of the mass transfer can easily be seen 
from the secularly changing period, its effect on the physical 
parameters of the system (e.g. masses and separation between 
the components) is below the accuracy of Ught curve modeUng. 
Therefore, we can assume that these values have been constant 
over the past fifty years covered by the observations. Hence, it 
should be legitimate to fold all the observations into a single 
phase diagram using the quadratic ephemeris (Eq. 4). We used 
the value of the period Pq, given in Table 3, for the reference 
epoch as the input to BINSYN. 

We combined the datasets (4), (7), (8), (9), and (11) from 
Table 1 into three groups corresponding to the Johnson's UBV 



Table 4. Initial and final values of the system parameters obtained with 
BINSYN. Errors of first four listed parameters were taken from Table 3 
and Eq. 6. 



parameter 


initial value 


converged value 


error 


Po [days] 


5.9044365 


- 


0.0000020 


Ml sin^ / [Mq] 


3.24 


- 


0.38 


a sin i [Rq] 


26.29 


- 


0.64 


g 


1.15 


- 


0.05 (est.) 


e 


0.0 


- 


- 


<•["] 


80.50 


80.51 


1.0 (est.) 


Hi 


8.77 




0. 10 (est.) 




4.00 


3.98 


0.02 (est.) 




1.0 






Fz 


1.0 






r,(pole) [K] 


13000 




1000 (est.) 


r2(pole) [K] 


5000 


5989 


200 (est.) 


M [jxM^y-'] 


3.60 


5.46 


0.50 (est.) 


Ra [Ro] 


9.2 


9.0 


1.0 (est.) 


Hv [Ro\ 


3.5 


2.9 


0.2 (est.) 



Table 5. Parameters of the UX Mon derived from the values listed in 
Table 4. 



parameter 


value 


Ml [Mo] 


3.38+0.40 


Mz [Mo] 


3.90+0.29 


i?,(pole) [/?o] 


3.49±0.05 


Riipole) [R^] 


9.80±0.03 


a [Re] 


26.65±0.65 


log gi (pole) 


3.88+ 0.01 


log g2 (pole) 


3.061±0.003 


Toivim) [K] 


6245+200 


rB(face) [K] 


7487+200 


/i 


0.39+0.03 


h 


0.44+0.01 




0.14+0.03 


/^(face) 


0.03+0.01 



filters and used the steepest descent method together with 
BINSYN to obtain the best fit of the system parameters. The 
solution was found by fitting the parameters of the system si- 
multaneously in all three passbands. 

Table 4 lists initial and final values of the system parameters. 
For the albedo and gravity brightening coefficients, we used the 
standard values for a radiative envelope of the primary and a 
convective envelope of the secondary. The initial values of pho- 
tospheric potentials (0| and Q2) and inclination, were chosen by 
trial and error. The same is true for the outer radius and height of 
the disk (Ra and Hy)- The iimer radius of the disk, Rb, was set 
to be equal to that of the star. 

Since we expect the secondary to be an evolved star that 
almost fills its Roche lobe, we assume synchronized rotation 
(Fi = !)• The situation is quite the opposite when we consider 
the primary. Because it might gain momentum from the accret- 
ing material, we expect it to rotate faster. However, initial trials 
showed that this parameter has no visible effect on the result- 
ing light curve. This is beacuse the primary is almost invisible 
through the accretion disk surrounding it. To reflect this, we set 
the formal value of this parameter Fi = 1, but its final value is 
completely uncertain. Our model does not offer the possibility to 
measure this parameter more accurately. 
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Fig. 8. Observed data and light curves calculated for the final set of Fig. 9. Three views of the binary system UX Mon with an optically 
parameters (Table 4). thick disk around the primary. 



The calculated and observed light curves in three passbands 
are shown in Fig. 8. In Fig. 9, we show three views of UX Mon 
at different orbital phases. The shape of the large secondary cor- 
responds to its Roche geometry close to the critical, Roche-lobe 
filling case. The primary star, which is mainly eclipsed by the 
disk, is outlined for enhanced clarity. 

We also investigated the uncertainty in the parameters de- 
scribing the geometry of the disk, Ra and Hy. The quality of the 
model fit for the various values of both parameters in the neigh- 
borhood of the best-fit solution is shown in Fig. 10. One can see 



that the uncertainty in the height of the disk Hy is much smaller 
than that of the outer disk radius Ra- Our estimates of the errors 
in both parameters, given in Table 4, are based on the mapping 
presented in Fig. 10. 

In Table 5, we give some important parameters of UX Mon 
derived from the values shown in Table 4. Note that, in both 
tables, certain errors are not indicated as estimated values, al- 
though they are calculated via a parameter whose error is esti- 
mated (mainly through connection with inclination, /, and mass 
ratio, q). Our estimate of the error in the inclination, i, is based 
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- - 0.00125 

- - 0.00120 

- - 0.00115 

10.00110 
0.00105 
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RaPoI 



Fig. 10. A map plot of rms vs. Ra and Hy in the vicinity of the light 
curve solution. 



on the work of Linnell et al. (2006), who also used the BINSYN 
software package in their analysis. We admit, however, that 
UX Mon is a more complicated system because of the circum- 
stellar matter around the primary. Consequently, the true uncer- 
tainty in the determination of i could even be several degrees. 
The estimate of the error in the mass ratio, q = 0.05, is subjec- 
tive, based on our experience with the sensitivity of KOREL in 
finding the mass ratio (Sect 3.2). We also admit that the param- 
eter mapping shown in Fig. 5 leaves us with some uncertainty 
regarding the error in q, which should not, however, be larger 
than ±0.1. The values of To are the temperatures of the disk and 
h, h, '^(rim), and Z/)(face) are luminosities relative to the system 
luminosity of the primary, secondary, rim of the disk, and face 
of the disk, respectively. Luminosities are calculated at phase 
(/> = 0.25 for a wavelength of A = 550 imi. 

4. Discussion 

Using the stellar parameters obtained from our light curve solu- 
tion, we can repeat the procedure of estimating the mass ratio, q, 
described at the beginning of section 3.2. By first using, JiCback) 
= 5742 K and log^2(back) = 2.83, we get ^, = (0.5 ± 0.25) km 
s"' and V2 sin; = 83.0 + 0.8 km s"'. The minimum obtained 
for this model is ~1.5, i.e. ~1.25 times smaller than the;^'^ ob- 
tained from the model described at the beginning of section 3.2. 

In Fig. 11, we show the observed and synthetic spectra as 
well as the difference between the two. The match is not perfect, 
but the deviations can be attributed to the disk that is still visible. 

We can again estimate the mass ratio, but this time using all 
values obtained in this paper: K2 = 104.6 km s~', i = 80.5", 
f = 5.9044365 days, and V2 sin / = 83.0 km s'K This time the 
estimated mass ratio turns out to be qest - 1 -22. This result is 
close to the value of q = 1.15 ± 0.05 that we found from the 
radial velocity analysis. 

By looking at the difl'erences between the measured and cal- 
culated RVs, we can use the Rossiter-McLaughlin effect to es- 
timate the size of the secondary. We assume that the phase at 
which the discrepancies start (0 =s 0.35) is the phase of the first 
contact of the secondary echpse. The phase interval between the 
first contact and the center of the secondary eclipse is then A<p « 
0. 15. The distance covered by the star during this phase interval 
is K2PA<p. On the other hand, we know from geometry of the 



1.3 




0.6 1 < < < < 1 1 1 1 — 

4900 5000 5100 5200 5300 5400 5500 5600 5700 

MA] 



Fig. 11. Observed (gray) and calculated spectra (black, superimposed 
over the observed spectrum) for 72 (back) = 5742 K and log g2 (back) = 
2.83. The difference spectrum is shifted upwards and shown in black. 
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orbital phase 

Fig. 12. Difference between the measured RVs and those calculated 
from Eq. 5. 



system, that this distance is equal to aiA^, where 02 is the dis- 
tance of the secondary to the center of mass. We can eliminate 
a2 by using 02 = R2/ sini2nA(p), assuming a circular orbit and 
ignoring the inclination. Finally we get 

K2P 

R2 = — sin(27rA0) ^ 9.9 Rq. (7) 
2n 

This value is close to the side radius of the disk that we found 
from our photometric solution /?2(side) = 10.3 Rq. From Fig. 13, 

we see that the eclipse starts at (f) x 0.365, which is a little later 
than we estimated from the RV measurements alone {(p ^ 0.35). 

5. Principal results and open questions 

The following principal results of our study lead to a substantial 
revision of the currently accepted model of the system: 

1. The orbital period of the system is secularly decreasing at a 
rapid rate of 0.26 s per year. 

2. The analysis of spectroscopic observations as well as the 
(plausible) assumption that the cooler, Roche-lobe filling 
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(|)=0.365 



Fig. 13. Phase of first contact during secondary eclipse found from pho- 
tometric solution. 



component has its rotation synchronized with the orbital rev- 
olution both lead to the conclusion that the cooler component 
is the more massive of the two. 

3. Reasonable basic physical properties of the components 
were obtained via modeling the light curves on the premise 
that there is an optically thick accretion disk around the hot- 
ter component, which is eclipsed in the primary minimum. 
We admit, however, that the remaining uncertainties in some 
of the derived model parameters are still quite large. 

4. The above findings identify UX Mon as a binary in a rare 
early stage of the mass exchange before the mass ratio rever- 
sal. 

There are, however, still a few open questions and possible 
alternative views that need to be solved and clarified using new 
dedicated observations and more sophisticated modeUng: 

- There is no really satisfactory and self-consistent explana- 
tion of the true cause of non-orbital light changes clearly 
seen, for instance, in Fig. 8. These variations were very sys- 
tematically investigated by Olson et al. (2009), who noted 
that the scatter in the light curve is largest in the phase inter- 
val 0.6-0.8 from the primary minimum, i.e. around the elon- 
gation with component 1 receding from us (see their Fig 2 or 
Fig. 8 here). These are the phases where the gas stream be- 
tween the components can be seen projected on to the disk 
of the primary star. They therefore tentatively concluded that 
the light variations are related to variations in the rate of 
the mass transfer, possibly induced by the magnetic activ- 
ity of the cool mass-losing star. There is, however, a very 
interesting plot in their Fig. 8 showing cycUc light variations 
at phases of totality during the primary minima with cycle 
lengths of some 800-1000 days. These variations are remi- 
niscent of similar cyclic variations found for other strongly 
interacting binaries like RX Cas (517 d; Kalv 1979), j3 Lyr 
(282 d; Harmanec et al. 1996) or AU Mon (417 d; Desmet 
et al. 2009). For all these three binaries, the light variations 
are probably strictly periodic and Desmet et al. (2009) argue 
that they are caused by variations in the circumbinary matter. 
All four systems, including UX Mon, are mass-exchanging 
binaries but only UX Mon is in the initial phase before the 
mass ratio reversal. 

- There is also another aspect worth considering: Harmanec 
et al. (1996) showed that the bulk of the Ha emission in 



yS Lyr does not originate from the optically thick disk but 
from bipolar jets, oriented perpendicular to the orbital plane 
and originating from the region where the gas encircUng 
the mass-gaining star hits the original gas stream fiowing 
from the Roche-lobe filling component (Bisikalo et al. 2000). 
Using only indirect arguments (mainly the phase offset of 
the RV curve of the Ha emission), Desmet et al. (2009) 
speculated about the possible presence of bipolar jets also 
for AU Mon. For the moment, this possibility cannot be ex- 
cluded even for UX Mon. 

- There are several possible ways to improve our understand- 
ing of UX Mon. First, spectrointerferometry of high res- 
olution could show what the character and geometry of 
the medium responsible for the observed Ha emission is. 
Continuing systematic spectral and photometric variations 
could reveal possible long-term changes in the Ha emission 
and its relation to the - already known - long-term light 
changes. A hydrodynamical modeling based on the binary 
parameters obtained in this study could also show whether 
our tentative model is compatible with the predictions from 
the theory. Finally, a more sophisticated modeling of the 
light curves including the effects of various possible cir- 
cumstellar structures, appears necessary, and not only for 
UX Mon. 

- Considering how crucial UX Mon might be to achieving a 
clearer understanding of the process of large-scale mass ex- 
change in this type of system, it seems obvious that the con- 
tinuation of its systematic investigation is worth the effort. 
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